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We present results for the mass spectrum and decay constants using non-perturbatively 0(a) improved 
Wilson fermions. Three values of and 30 different quark masses are used to obtain the chiral and 
continuum limits. Special emphasis will be given to the question of taking the chiral limit and the existence 
of non-analytic behavior predicted by quenched chiral perturbation theory (q\PT). 



1. INTRODUCTION 

The chiral and continuum extrapolations of 
lattice results in quenched QCD have become 
an important issue since more accurate data 
for small quark masses and results for different 
lattice spacings have become available. Our 
simulations have been done for three different 
values of the gauge coupling, (3 = 6.0, 6.2 and 

6.4 (a « 0.09 — 0.05 fm), and on lattices of size 

1.5 — 2.3 fm. Our propagators have been gen- 
erated for 7-12 different values of the hopping 
parameter n with the ratio of pseudoscalar and 
vector meson mass mps/my in the range of 
0.41-0.98. 

2. CHIRAL EXTRAPOLATION 

Based on small quark mass expansions, chi- 
ral perturbation theory can give information 
about how the extrapolation to the chiral limit 
has to be done. For quenched QCD, qxPT 
predicts non-analytic terms, e.g. for the pseu- 
doscalar mass it gives to one-loop Q 



For full QCD, the parameter 5 should be zero. 

In order to determine S we fitted our data 
using the ansatz 

am PS /rhf = Cp (arh^)-^s + A x # amf 
+ A 2>0 (am* ) 2 ] / [1 + B hP afhf] , (2) 

where = Z M (go) ^^™^ m q is the 

renormalization group invariant, and rh q the 



2.4 



2.0 



1.6 



E 

ra 1.2 



• |3=6.0 

♦ P=6.2 
A P=6.4 



TOps/rhq = Corhq 



Clfh q 



0{rh 2 q ). 



(1) 



*Talk given by D. Pleiter at Lattice '99, Pisa, Italy. 
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Figure 1. The ratio rn-ps/^q as a function of 
rhq 1 . The solid line is a fit using eq. (||). 
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Figure 2. Vector meson mass as a function of 
the pseudoscalar mass. The solid line is a fit 
using eq. (Q) with N = 4. The dashed line is a 
fit using a cubic polynomial with amps < 0.6. 

bare Ward identity quark mass. The non- 
perturbative renormalization constant and the 
perturbative improvement coefficients 6^ and 
bp have been computed in ||. According to 
this ansatz rnp s /m^ should be linear up to 
heavier quark masses. Our results are shown 
in Fig. 0. From this fit we get S « 0.14(2). 

For the vector meson mass my and the nu- 
cleon mass tun, qx^T predicts |3|Q] 

amv,N — aMy,N + C1/2 amps (3) 
+d (amps) 2 + O ((amps) 3 ) 

with C'1/2 ~ —1.5 for my and C1/2 = 
-(3tt/2) (D - 3F) 2 (5 « -0.24 for the nucleon 
(D and F taken from f§). For full QCD, C 1/2 
should be zero. 

To reduce the number of free parameters in 
our fits, we make the assumption that the co- 
efficients C1/2, Ci, ... do not show any scaling 
violations, i.e. do not depend on f3. To take 
the different scales into account we replace in 
eq. (||) a — > ap = a0 = Q.oSp and use the scale 
parameter sp as an additional fit variable. We 
fitted our data using a Pade-like ansatz, 



apM V:N + Y*=i A i {apmps) 1 



apm v ,N 



1 + Bn-i (apmps) 



N-l 



(4) 



which is again chosen to give the correct heavy 
quark limit, and a polynomial ansatz, which 
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Figure 3. Fits as in Fig. || but with C1/2 = 0. 
The dotted line comes from a fit using eq. (g). 

corresponds to Eq. [| with fljv-i = 0. For 
the Pade-like ansatz the results for My,N are 
close to constant for 0.1 < Bjv-i < 1-5, we 
therefore use Bm-i = 1. For my we find 
C1/2 = -0.6(1) and -0.7(3) for the Pade- 
like and the polynomial ansatz, respectively. 
The results are plotted in Fig. |^. For mjv we 
find Ci/2 = 0.8(8) and 0.4(5). Our results for 
sp deviate at most 2% from the correspond- 
ing values obtained from the force scale r || , 
which can be explained by scaling violations. 

The fits shown in Fig. ^| can be compared 
with fits shown in Fig. ||, where C1/2 has been 
set to zero. Also shown are the results for a 
phenomenological ansatz 

3 

E 

i=2 



{a m v , N ) 2 = (apM VjN ) 



A^apmpsfib) 



which has less free parameters. 

3. CONTINUUM LIMIT 

For improved fermions we expect scaling vi- 
olations to be of 0(a 2 ). We therefore extrap- 
olate linearly in (a/ro) 2 . To avoid artifacts of 
the quenched approximation, we use the phe- 
nomenological ansatz (amv^jv) 2 = (aMv,N) 2 + 
A2,[3(amps) 2 + A3 t p(amps) 3 for the extrapo- 
lation to the chiral limit, while we use a Pade- 
like ansatz similar to eq. (|J) with Ci/ 2 =0 to 
interpolate our data or to extrapolate towards 
the charm quark mass region. 
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Figure 4. Continuum extrapolation of the vec- 
tor meson masses. The stars show the experi- 
mental value. 
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Figure 5. The same as in Fig. |4] for the pseu- 
doscalar decay constants f„ and fx- 

In Fig. ^ we show our results for the vector 
meson masses. We used the experimental val- 
ues for the 7T, K, and D mass to fix n Ut d 1 k s , 
and K c h, respectively. Since we use results for 
degenerate quark masses, we have to restrict 
ourselves to the light masses in the baryon sec- 
tor. The results are shown in Fig. ^. 

To one-loop q\PT, the pseudoscalar de- 
cay constants with degenerate quark masses 
should not suffer from quenched artifacts. We 
again used the phenomcnological ansatz to ob- 
tain the chiral limit. The continuum extrapo- 
lation of / w and Jk is shown in Fig. ||. 

4. CONCLUSIONS 

We found evidence for non-analytic behav- 
ior of m 2 PS /m q and my. The results for 
to TV do not seem to confirm the predictions 
of q\PT, but within errors we cannot con- 
clude Ci/2 > 0. Our result for the parameter 
S agrees with the results of the CP-PACS col- 
laboration and is slightly larger than the 
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Figure 6. The same as in Fig. [| for spin-| 
octet baryons. 

number found by the FNAL group Q . For my 
we found Cj/2 to be much larger than in 0. 

If we use the phenomenological ansatz for 
the extrapolation to the chiral limit the ratio 
TO/v/TOp seems to become larger than the ex- 
perimental value. This is probably a quench- 
ing effect. For all masses and decay constants 
presented here discretization errors are 0(a 2 ). 
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